ASRC Aerospace Corporation 


P.O.Box 21087 

Kennedy Space Center, Florida 32815-0087 

SBA Section 8(a) Company 
Alaska Native Corporation 


NIAC CP 04-01 
Phase I 

Advanced Aeronautical/Space Concept Studies 
Analysis of a Lunar Base Electrostatic Radiation Shield Concept 


Mid-Term Phase I Report 


Charles R. Buhler, Principal Investigator 
(321) 867-4861 
Fax: (321) 867-4489 
Charles.Buhler-1 @ksc.nasa.gov 


Leon Wichmann, Contracting Officer 
(321) 867-1504 
Fax: (321)867-1087 
Leon. Wichmann- 1 @ksc.nasa. gov 


December 15, 2004 



Introduction 

Space weather can be defined as the total ensemble of radiation in space, as well as on the 
surface of moons and asteroids. It consists of electromagnetic, charged-particle, and neutral 
particle radiation. 

The fundamental goal behind this NIAC Phase I research is to investigate methods of generating 
a static electric-field potential ®(x,y, z ) in the volume above and around a “safe” or protected 
area on the lunar surface so that trajectories of harmful charged particle radiation are modified 
(deflected or reflected), thus creating a shadow over that region. Since the charged particles are 
not neutralized but merely redirected, there will be areas outside of the shadowed protected 
region that will have a higher flux concentration of radiation. 

One of the fundamental limitations of the static electric (electrostatic)-field approach to radiation 
shielding is that complete shadowing is accomplished only by complete reflection, which can 
only occur for shield voltages greater than or equal to the kinetic energy (in electron volts) of the 
incoming charged particles. Just as habitats on Earth are protected from severe weather events 
and conditions, such as extreme temperatures, high winds, and UV radiation, using multiple 
methods of shielding protection from severe space weather will undoubtedly require multiple 
strategies. The electrostatic shield concept may be one of many methods employed to protect 
astronaut habitats on the lunar surface from some of the harmful effects of space weather. 

Two Design Concepts: Electrostatic Spheres and Electrostatic Screens 

The goal of this Phase I award is to investigate electrostatic field configurations and methods of 
generating those fields that will reduce the intensity of charged-particle radiation on the lunar 
surface in a protected volume. A major design challenge in this approach is the requirement to 
shield both negative and positive particles. Complete shadowing of radiation (total reflection) is 
the ideal goal, but partial reduction in radiation flux may also be useful if the electrostatic shield 
is one component of a multistage radiation protection system. 

Two parallel design activities have been pursued concurrently during the first half of Phase I. 
ELECTROSTATIC SPHERES 

The motivation behind this approach is that generation of an arbitrarily specified electrostatic 
field can be approximated by a system of conducting spheres of specific voltages and diameters. 
A distinct advantage of computing the 3-D electrostatic-field potential from a system of 
conducting spheres is that direct analytical methods can be used, thus significantly reducing the 
computation complexity. 

A close analogy exists to this design methodology in the field of audio engineering: a parametric 
equalizer is used to approximate a user-defined frequency response H(co) by summing the 
individual responses of second-order band pass filters of variable bandwidth, center frequency, 
and amplitude. If the electrostatic-field potential Q>(x,y, z ) is plotted along a line defined by the 
parameter 5 , ® (As + a, Bs + b, Cs + c ) , the resulting plot of ®(s) versus s would look similar to a 

plot of H(cb) versus co. Well-known design techniques of finding optimized parametric filter 
parameters exist in order to approximate a user-defined frequency response. 

Similar mathematical optimization techniques could be applied to finding the best-fit 
electrostatic-sphere parameters that match a desired 3-D potential contour in the radiation shield 
design problem. This may be an area of future work, but in this present activity, most of the 
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effort was devoted to generating and testing software to calculate the potential field and perform 
Monte Carlo simulation of charged-particle trajectories, with and without the powered spheres. 
The design configurations of Phase I have been based on empirical and best-guess methods of 
specifying the sphere parameters. Details of this method are discussed in Appendix A. 

ELECTROSTATIC SCREENS 

A finite-element approach prevents us from falling into the trap of advocating a design solution 
simply because it is easy to calculate. Field Precision (New Mexico), our Phase I collaborator, 
has developed a suite of software modules to compute electric fields, magnetic fields, and 
charged-particle trajectory analysis, based on finite-element modeling. Field Precision is using 
% its software to investigate a second method of radiation shielding based on electrostatic screens. 
Details of this method are discussed in Appendix B. 

Electrostatic Shield Design Problems and Constraints 

The electrostatic shield is composed of conducting surfaces (or electrodes), connected by 
insulating components, and mechanical support structures. The shield may be composed of 
spheres, screens, combinations of the two, or something else. The goal is to find the best overall 
solution for creating an ideal static electric-field potential, based on electrical, mechanical, and 
structural engineering, as well as economic constraints. Two parts to this problem are difficult to 
separate: (1) define the ideal electrostatic field for radiation shielding, and (2) find a way to 
generate it. 

Assuming that generators can be fabricated that produce the necessary voltages to power the 
electrostatic radiation shield, the fundamental design constraints of the system that must be 
considered are: 

• Insulating structures, such as supports, cables, and poles, must be limited to an electric- 
field strength exposure of no more than some maximum value, Evb in order to avoid 
vacuum breakdown. A conservative value for this maximum might be 2.5 [MV m" 1 ]. 

• Sunlight-generated electron photoemission from the lunar surface must be blocked from 
interfering with the shield. 

• Lunar dust must be blocked from contacting with the shield’s charged surfaces. 

• Charged-particle radiation from the horizon will most likely not be stopped by any 
practical shield configuration — a wall of lunar sand bags around the base may be a good 
solution. 

• Attractive and repulsive coulomb forces between spheres, electrodes, and all other high- 
voltage structures must be balanced by appropriate mechanical support structures. 

• Size and mass of deployable structures must be limited to what can be assembled by 
astronauts on the lunar surface. 


Shield Efficiency and Total Radiation Dosage Reduction 

For the purposes of radiation shield analysis, it is useful to construct a model quantifying the 
dosage rate of harmful radiation received by an unprotected astronaut on the lunar surface (or in 
deep space). This can be empirically modeled as: 
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A>(0 = \p(E,t) /3(E) dE 

o 

where p(E,t) is the number density of charged particles that intersect a critical surface, or 
protected region. This critical surface could correspond to the entire surface of an astronaut’s 
body or some other arbitrary surface surrounding the astronaut’s habitat. The double integral of 
p(E,t ) over all energies E and time t, gives the total number of particles n 0 , a dimensionless 
quantity, intersecting the critical surface. 

In order to estimate the biological damage from a particle of energy E, a damage coefficient 
/ 3 (E) can be used to approximate this effect. The many simplifying assumptions called upon 
here include ignoring the particle angle of incidence, atomic mass of the particle and 
composition, charge of the particle, and specific region of the body impacted by the radiation. 
The only parameter used in this simplified model is the particle kinetic energy, E. /3(E) is the 
relative amount of biological radiation damage as a function of energy. (3(E) may be expressed 
as a polynomial series, a power-law, or some other approximation function. 

A shielding efficiency e(E ) can be defined (see Figure 10 of Appendix B) so that the total 
dosage rate of harmful radiation is reduced as described by the following integral: 

co 

D(t)m \{\-e(Ej\p(E,t)/3(E)dE (2) 

0 


For a shield of perfect efficiency (i.e., s(E) - 1) for all E, D(t) = 0, according to Equation (2). 
Realistically, the fraction of accumulated harmful radiation dosage can be expressed as the ratio 
of the time integrals of Equations (1) and (2): 


\D(t)dt 

f ° ' j D„m 


£ D — 1 f D 


(3) 


where we define as the dosage attenuation efficiency. In the case of a specific flux of 
radiation intersecting a critical surface, the specific efficiency £s can be defined by calculating 
the total number of intersecting particles with (ns) and without (no) the shield: 

fjtl-<r(£)KE,OaK<ft 
5 \\p(E,t)dEdt 


The specific efficiency £s of Equation (4) will generally show a more optimistic quantity when 
compared to the fraction of accumulated radiation dosage £d since (3(E) would be expected to 
increase with energy. Only in the special (and unrealistic) case of (3 (E) = 1 , would % D =% S , but 
normally, £ 0 <£ s . 
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Future Work 


Topics for future work are listed below. Some of these topics may be appropriate as part of a 
possible Phase II, whereas other topics listed will be completed during Phase I. 

• Future simulations will consider both a ground shield for protection against lunar surface 
dust and photoemission electrons, as well as an outer electron shield to deflect incoming 
electron radiation from space. It may be possible to combine the negative shield and the 
ground shield into one structure. 

• A combination of screens (concept #2 - Appendix B) and spheres (concept #1 - 
Appendix A) may be advantages. For example, the negative ground shield would be 
composed of an electrostatic screen, whereas the high-voltage shield would be a set of 
+100 MV spheres on 40-meter poles. 

• Electrostatic rings or “donuts” have been suggested. These could be configured numerous 
ways, such as concentrically, forming a cone over the protected area. 

• The big challenge with the electrostatic-shield approach is the incredibly large voltages 
(an order of magnitude larger than anything obtained to date) needed to deflect ionized 
space radiation. Future work will investigate a concept to reduce the “base” field strength 
by using small transverse magnetic steering fields. The electrostatic shield works 
somewhat like a cathode ray tube (CRT), where the anode is analogous to the shield and 
where the cathode generates the radiation source. In the case of the electrostatic shield, 
the anode is meant to deflect oncoming charges particles. A CRT incorporates low- 
current magnetic steering fields. The goal of this activity is to investigate the possibility 
of applying the magnetic steering field concept to the electrostatic shield in order to 
reduce the large electric fields needed. 
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Appendix A: Electrostatic Sphere Concept 

The following report describes a FORTRAN simulation model, originally developed at Kennedy 
Space Center, NASA Applied Physics Lab (2002 to 2003), for modeling charged-particle 
radiation scattering from electrostatic spheres. This report describes the mathematical details, 
based primarily on modeling the dynamics of a charged particle in an electric field. The particle 
velocity is determined using relativistic mechanics, from user specified kinetic energy 
distributions. The primary modifications made to that software for the NIAC Phase I project are 
as follows: 

• Incorporated a zero-potential ground plane, representing the lunar surface, using image 
charges below the surface. 

• Incorporated a hemisphere isotropic radiation trajectory distribution, representing the 
lunar sky. 

• Incorporated two particle Monte Carlo energy distributions, for two independent 
particles of arbitrary mass and charge. This allows simultaneous study of electron and 
positive-ion radiation. 



The above concept artwork depicts an electrostatic-sphere-based radiation shield. In this 
example, the lower ground shield is at a medium voltage potential (approximately -100 kV), 
prevents electrons from space from reaching the surface, and prevents electrons from the surface 
from reaching the upper high-voltage spheres. The walls, composed of bagged lunar regolith, 
stop charged particles from the horizon. Note that this figure is for conceptual purposes only and 
does not necessarily scale properly. For example, if the high-voltage spheres were at +100 MV, 
then the poles would need to be 40 meters high to satisfy the constraint that maximum fields are 
kept below the vacuum breakdown value of 2.5 MV m' 1 . 
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Figures 3 through 6, following the Appendix A report, are simulations of a complex sphere 
configuration consisting of seven upper negatively charged spheres and three smaller positively 
charged spheres. The protected region is arbitrarily defined as a hemisphere of 4-meter radius at 
x = 0, y = 0, z = 0 (on the lunar surface). The software generates 10,000 particles, half of which 
are 50-MeV protons and half are 5-MeV electrons. The first of each set of figures show the 
simulation results with zero shield voltage (i.e., no external electric fields). This provides a 
comparison to powered shield shown as the second figure of each set. For this particular 
simulation example, the number of particles that intersect the protected region in the unpowered 
case is N off = 323 and N on - 20 in the powered case. Therefore, specific shield efficiency & from 
Equation (4) can be written as 

&=1 (A-l) 

N 

^ Off 

In the case of this simulation, & - 0.94 or 94 percent. 
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THE ELECTRIC Field due to a system of N point charges at a field point r is given 
by: 


E(r) = 


1 

AttSo 



( 1 ) 


where r, is the location of the ith point charge q . If the /th point charge is implemented 
as a sphere with radius and a uniform charge distribution at potential V j , Equation (1) 
can be rewritten as: 


E(r) = 2M 


/=! 


r-r 


r-r. 


( 2 ) 


THE TRAJECTORY of a charged particle of rest mass m 0 and charge Q in the electric 
field given by Equation (1) is determined by its equation of motion: 


0E(r) = p = — (y m 0 \) 

at 


= ym { 


f , V • V '' 
- V 


v + y 


where y = (l - v 2 / c 2 ) 


( 3 ) 


The acceleration of the particle, a = v , of the particle is calculated by re-writing Equation 
(3) as: 


where, 
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Solving for a in Equation (4), 


where, 


and, 


a = — C _ 1 E(r) 
ym o 
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B — (^13^22^31 C^CjjCji CjjCjjCjj + C| 1 C 23 C 32 + C|2^2I^33 ^ 1 1^22^33 ) Y 

A x — B x < -'22 ( '33'^.t _ C l3 C 32^y ^12^33 ^y < 'I3 < '22 ~ ^12^23 ) Q 
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A z = (c 22 c 3 jE x — c 2] c 32 E x —c 12 c 3t E y +c n c 32 E y + c ]2 c 2] E z — c u c 22 E z }Q 


( 6 ) 


(7) 

(8a) 

(8b) 

(8c) 


Based on a Taylor series expansion about time point k, a set of difference equations for 
position and velocity can be expressed as: 


v* + i * v, + \ k At 
*v t + a, At 


(9a) 


r * +1 «i' i +f i A/ + if,At 2 
*r, + v jt A? + 1 a, At 2 


(9b) 


In the non-relativistic case, the acceleration in Equations (9) reduces to a k =(Q/ m 0 )E k . 
In the general case (non-relativistic and relativistic), a k is computed from Equation (6) 
using r k , v k , and E(r k ) . 
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INITIAL State Vector of the particle is modeled by setting the initial particle 
position r 0 at a constant far field distance R' 0 from the vehicle (dotted spherical surface 

in Figure 1), using random parameters, </>' and 0’ . The particle’s initial velocity v 0 is 
composed of an initial speed v 0 and a direction. The direction of v 0 is given by the two 
random variables, as shown in Figure 1, angles <p and 6 . The vehicle is fixed at the 
origin of the coordinate system and therefore the particle’s initial state vector and 
trajectory is expressed relative to the vehicle’s coordinate system. The software is 
configured to reject all initial states whose initial trajectory vector does not intersect with 
the inner sphere of radius R 0 . This is determined by performing the following 

comparison: r 0 + .R' v 0 |v 0 | - ' <R 0 . If false, the trajectory is rejected and a new one is 
calculated. 
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The particle’s initial velocity vector is given by the following: 



/ \ 





a 

v o = 

V , 

= v o 

p 


l V J 




where, a 2 + j3 2 +y 2 = 1 , and: 


(10a) 


a = sin^cos^ /?ssini9sin^ yscos 0 (10b) 

The particle’s initial velocity v 0 does not change until the particle encounters an external 

force. For the purpose of simulation, it is convenient to assume that the force from the 
vehicle’s electrostatic shield has no influence on the particle’s trajectory until it is within 
a radius R' 0 from the vehicle. The initial position vector r 0 is (see Figure 1): 


r o ~ *o 


^sin O' cos <f>'^ 
sin O' sin <j>' 


( 11 ) 


^ cos O' J 

An option has been added in version 1.0: the semi-sphere radiation flag which limits the 
radiation trajectories to the upper hemisphere, for the purpose of lunar base version. 


PROBABILITY DISTRIBUTIONS of the individual model parameters are based on the 
uniform distribution P u equal to one in the interval from 0 to +1, and zero elsewhere; a 

sinusoidal distribution P s (u ) = sinw in the interval from 0 to +1, and zero elsewhere; and 
the normal distribution P N of unit variance and zero mean, defined by: 


PM = 



(12a) 


Table 1 describes the individual model parameters and their corresponding distribution 
functions. R' 0 is the distance to the initial position of the charged particle, while p 0 is 

the critical radius of influence. As shown in Figure 1, the angles 0 and O' describing the 
incoming particle velocity, range from 0 to n , whereas the angles <f> and f range from 0 

to 2n . Note that R 0 < R' 0 and p 0 <R' 0 , however, either p 0 < R 0 or p 0 > R 0 is valid. 


The speed v 0 of the particle is determined from the kinetic energy, T : 


cj_ l + 2£ 

i+4 


4 - 


To 


(13) 


where m 0 is the rest mass of the particle. The energy of the particle is detennined by the 
standard deviation <j e and mean T 0 . 
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Table 1. Model Parameters and Associated Distributions. 


/ 


Model 

Parameter 

Associated 

Constant 

Type of 
Distribution 

Distribution 

Formula 

Po 

Po 

constant 

- 

R o 

R o 

constant 

- 

K 

K 

constant 

- 

9 

K 

sinusoidal 

7TP S 

9' 

n 

sinusoidal 

nP s 

<t> 

In 

uniform 

2nP u 

</>' 

2 n 

uniform 

InPy 

T 

T 0 

normal 1 

+ T 0 I 


' The charged particle energy distribution will be replaced with a more 
meaningful distribution function in future versions of the model. 


If the energy of the particle is expressed in electron volts [eV], the speed v 0 of the 
particle from Equation (13) is modified as: 


where the mass and velocity are expressed in mks units, and q is the fundamental unit of 
charge * 1.6 10" 19 [C], 


Plotting of Particle Trajectories: in order to effectively plot the particle 
trajectory and intersections with the critical radius sphere defined by p 0 , three 

projections of the 3D volume particle paths are generated as xyPlot, xzPlot, and yzPlot. 
These plots are not as easily interpreted as the previous plots where the trajectories 
intersected with a disk. However, after some thought, it should be easy to understand 
what is occurring. 

Previous versions of ESM used p Q as the plot scale factor, such that the plot extent both 
in the horizontal and vertical directions, is 2 p 0 . Starting with this version ESM (v5.1), 
and all future versions, the plot will be scaled by a separate parameter, p p , such that the 
plot extent will be 2 p p in both dimensions. 
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IMAGE Charge CORRECTION applied to the system of conducting spheres, results in 
a configuration as depicted in Figure 2. 


z • zeroth order charge, 



Figure 2. Image Charge Treatment of a System of Conducting Spheres. 

The electric field from Equation (1) now becomes a double summation over all N 
spheres, where there are N 2 total charges, N zeroth order charges and N(N - 1) first 
order image charges: 


1 N 

E(r) = Y 

4^ 0 tr 
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(15) 
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where the image charge parameters are given by the following relations: 
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The potential on the surface of sphere k can be expressed as: 
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(17) 


j*> 


J 


where e is any unit vector. Since all charges are unknowns, the goal is to eliminate all q’s 
from Equation (15). 
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Equation (18a) can also be expressed in matrix form as: 

V = Cq 

and solving for q: 

q = C’ 1 • V 


(18a) 


(18b) 


(19a) 


(19b) 


< 
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where V is a vector of length N of sphere voltages V k , q is a vector of length N of zeroth 
order charges , and C is an Nx N array as described by Equation (18b). 

Since e is a random unit vector, C according to Equation (18b) will not have a unique 
value. A better estimate of C is obtained by averaging Equation (18b) over many random 
unit vectors, e/: 


Q,= 


1 1 


M 


M£t 


N 

I 


R j 


\ R k*i + r t -*i| P R k *i + r, - r j, r j 




( 20 ) 


The weakness of this method may be related to the ambiguous solution which is 
dependent on the specific choices of e/ in Equation (20). As the number of values 
averaged over increase (i.e., as M increase), the solution given by Equation (19b) may 
converge to a unique value. 
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Figure 3a. Simulation Run of Unpowered Lunar Electrostatic Shield Model (LESM vl.2)~ User 

Interface and Sphere Configuration File. 
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Figure 3b. Simulation Run of Powered Lunar Electrostatic Shield Model (LESM vl.2). 
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Figure 4a. Simulation Run of Lunar Electrostatic Shield Model (LESM vl.2) - Red dots are intersection 
of electrons and blue dots are intersection of protons with lunar surface (z = 0). Gray circles are x-y 

projections of unpowered electrostatic spheres. 
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Figure 4b. Simulation Run of Lunar Electrostatic Shield Model (LESM vl.2) - Red dots are intersection 
of electrons and blue dots are intersection of protons with lunar surface (z = 0). Purple and cyan circles 

are x-y projections of powered electrostatic spheres. 
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Figure 5 a. Simulation Run of Lunar Electrostatic Shield Model (LESM vl.2) - Red dots are intersection 
of electrons and blue dots are intersection of protons with a 4 [m] radius sphere (protected area) 
centered at x = 0, y - 0, z = 0. Gray circles are x-y projections of unpowered electrostatic spheres. 

Yellow-gray trails are particle trajectory paths. 



Figure 5b. Simulation Run of Lunar Electrostatic Shield Model (LESM vl.2) - Red dots are intersection 
of electrons and blue dots are intersection of protons with a 4 [m] radius sphere (protected area) 
centered at x = 0, y = 0, z = 0. Purple and cyan circles are x-y projections of powered electrostatic 
spheres. Yellow-gray trails are particle trajectory paths. 
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Figure 6a. Simulation Run of Lunar Electrostatic Shield Model (LESM vl.2) - Red dots are intersection 
of electrons and blue dots are intersection of protons with a 4 [m] radius sphere (protected area) 
centered at x = 0, y = 0, z = 0. Gray circles are y-z projections of unpowered electrostatic spheres. 

Yellow-gray trails are particle trajectory paths. 
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Figure 6b. Simulation Run of Lunar Electrostatic Shield Model (LESM vl.2) - Red dots are intersection 
of electrons and blue dots are intersection of protons with a 4 [m] radius sphere (protected area) 
centered at x = 0, y = 0, z = 0. Purple and cyan circles are y-z projections of powered electrostatic 
spheres. Yellow-gray trails are particle trajectory paths. 
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Appendix B: Electrostatic Screen Concept 

The following report is a detailed analysis of an electrostatic radiation shield design, done 
independently by Field Precision (New Mexico). Some of the highlights and conclusions of this 
report are: 

• The theoretical best shielding efficiency (where the maximum is 1) goes as: 

£(£) = V<?<V£ (B-l) 

where ®o is the maximum shield voltage, q is the charge on the particle, and E is the 
kinetic energy of the particle in electron volts. 

• Field Precision software simulations support the conclusion that electrostatic screens 
provide a more efficient shielding system, one that approaches the ideal efficiency given 
by Equation (B- 1 ) 

• A negative ground shield is essential to stopping electron photoemission. The ground 
shield may also be effective at solving the lunar dust contamination problem. 



The above artwork figure illustrates the electrostatic screen shield concept. Again, this figure 
is not to scale and comments similar to those in Appendix A apply (i.e., maximum fields 
must be kept below the vacuum breakdown value so that the poles would need to be 40 
meters high in order to achieve a useable shielding voltage of ®o = 100 MV. Also, the mesh 
density of the screens would be much higher than that illustrated, and the area extent of 
screen coverage would need to be greater in order to shield the working area shown in the 
figure. As with the previous sphere concept described in Appendix A, the screen 
configuration would not be able to stop low-elevation-angle radiation from coming in from 
the horizon; thus a wall would be needed for complete protection. 
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1 Shield system concept 

It is important to delineate a clear set of goals for and constraints on an elec- 
trostatic shielding system for a moon base before proceeding to an extensive 
program of three-dimensional simulations. To begin, I shall list several basic 
facts. The surface of moon is approximately an equipotential plane. In the 
following discussions, I shall take the surface potential as a reference (ie., 
ground potential equal to 0.0 V). Clearly, personnel and equipment on the 
moon surface must be approximately at ground potential. The potential of 
the infinite space above is also approximately 0 = 0.0 V. Any difference in 
the potential of the moon surface and space would result in an exchange of 
low-energy electrons that would restore the balance. The issue of concern is 
the flux of energetic ions that arrive from space as a result of solar proton 
events or as part of the galactive cosmic radiation. The spatial distribution 
of incident ions is approximately isotropic. The role of an electrostatic shield 
is to reduce the energy-flux of ions that strikes a protected area on the moon 
surface. No matter what the electrode geometry, the basic function of the 
shield is to generate a region of high positive potential above the protected 
area that will reflect a high proportion of incident ions. Note that the shield 
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cannot change the average energy of non-reflected ions that are able to reach 
the surface. These ions start in space at the reference potential and strike 
objects at the same potential. I discussed a second constraint in my first 
report - electric field components transverse to the direction of ion motion 
will not change the average flux of unreflected ions if the incident distribu- 
tion is isotropic. For a given applied voltage, the implication is that the 
most effective shield is one that maintains a uniform positive potential over 
the protected area. Therefore, in this report I consider sheet or mesh elec- 
trodes. Multiple electrode arrays would only yield an advantage if they had 
significantly smaller weight. 

Figure 1 shows the preliminary concept for a moon base shield presented 
in our proposal. This approach has several problems: 

• The electrodes over the ground plane are spaced too far apart compared 
to their height. The magnitude of the potential barrier at the center of 
the array will be only a fraction of the applied voltage on the spheres. 

• The conceptual design of the high-voltage system is optimistic. The 
figure shows spheres on thin poles with a nearby small box to create 
the voltage. In actuality, each sphere would rest on a Van De Graaf 
generator operating at 20 MV or higher, a piece of equipment that 
would dwarf any of the structures on the ground. Disregarding the 
issue of whether a multi- MV column can operate in vacuum, the use of 
multiple generators would be wasteful. 

• Intense electric fields on the lunar surface would polarize dust particles 
which would be attracted to regions of high electric field gradient. 

• Personnel and equipment are not protected against strong fields, so 
there is a possibility of a damage or electrocution in the event of a 
breakdown. 

The next section describes an alternate geometry that resolves the problems. 

2 Improved shield system 

Assuming that high-voltage engineering problems can be solved, Figure 2 
shows an improved version of an electrostatic shield. The system utilizes a 
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Figure 1: Conceptual electrostatic shield for a moon base from the proposal 
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Figure 2: Diagrammatic view - electrostatic shield system for a moon base 

single, centralized voltage source (possibly in the form of a black monolith 
to blend in). For an applied voltage of 20 MV and a conservative gradient 
of 2.5 MV/m for the vacuum insulators, the generator would rise about 8 
m above the ground shield. The high-voltage shield consists of a thin wire 
netting supported on insulating poles. Voltage grading structures on the 
outer edge would reduce the chance of breakdown. It may be necessary to 
cover the central portion of the netting with an opaque material to reduce 
photoelectron emission on the insulating column. 

A larger second net with a negative applied potential would be suspended 
about 3 m from the surface of the moon. The ground shield serves several 
important functions: 

• It reduces the generator load caused by photoelectric emission from the 
ground. 

• It acts as a Faraday cage to protect personnel and equipment from 
breakdowns and photoemissive charging. 

• The shield reduces the field at the moon surface to prevent attraction 
of dust particles. 
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Regarding the first item, the significance of photoelectron loading is reviewed 
in Sect. 3. Section 4 addresses required voltage levels on the ground shield 
to ensure a net field at the surface that repels electrons. 

3 Electron photoemission from the lunar surface 

A moon base would be exposed to sunlight half the time. The ultraviolet 
portion of the solar spectrum can generate photoelectrons on the lunar sur- 
face as well as on equipment and personnel. Ordinarily the emitting objects 
would be left with a small positive potential (a few volts) that would attract 
the low-energy electrons to preserve charge neutrality. Without the ground 
shield of Fig. 2, the HV electrode would generate a strong attractive field 
on the lunar surface. In this case, all photo-electrons would travel to the 
HV shield, loading the voltage generator. Even though the current density 
may be small, the available emission area is very large. In the configuration 
discussed in Sect. 5, the area exceeds 3000 m 2 . A related problem is the 
selective charging of personnel and equipment if their work function differs 
significantly from the lunar surface. The emission rate from metal equip- 
ment could be substantially higher than lunar dirt, resulting in large voltage 
differences and possibly arcs. 

We can make a rough estimate of the available current density to gauge 
the severity of the problem. First, consider the available photons. Figure 3 
plots the integrated solar power flux as a function of photon energy (hi/) in 
eV. The data were adapted from Ref. [1]. The plotted quantity is the integral 
of power flux over all higher energies. The second piece of required informa- 
tion is the electron yield as a function of hi/. Here, the term yield refers to 
the number of emitted electrons per incident photon. There is considerable 
variation between materials. Fig. 4 shows data abstracted from Ref. [2] for 
copper and silicon surfaces, while Fig. 5 shows theoretical predictions for 
emission from small dust particles from Ref. [3]. It is difficult to make an 
accurate calculation because the solar spectrum drops rapidly in the energy 
range of rising emission coefficient. To make a conservative estimate, we take 
an average yield of 10 -3 for the portion of the photon spectrum above 6 eV. 
The photon flux above 6 eV is approximately equal to 


F p = (0.02 x 10 3 )/ (6) (1.6 x 10 19 ) = 1.25 x 10 14 (photons/s/cm 2 ). (1) 
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copper with clean surface [18]. 



Figure 25.13 PE quantum yield of n- and p-type 
silicon with different doping {x/E g ~ 4) [16]. 


Figure 4: Photoelectric yield as a function of photon energy for copper and 
silicon 

The available electron current density is 

j e = (10 _3 )(F P )(1.6 x 10-19) = 2.0 x 10 -8 (A/cm 2 ). (2) 

Multiplying by an area of 3.0 x 10 7 cm 2 gives a leakage current of 0.6 A. 

To maintain voltage, the shield would require a large input power (ex- 
ceeding 12 MW). Clearly a practical system must include the ground shield 
of Fig. 2. Photo-emission from the shield can be minimized by using a mesh 
with high transparency and including a coating with high work function. One 
remaining problem is the attraction of electrons from the space environment 
above the HV shield. I will collect data to make estimates in a following re- 
port. If there is significant electron current from space, it may be necessary 
to include an additional ground shield above the HV shield. 


4 Ground shield design considerations 

I envision the ground shield as a high-transparency mesh formed from a 
square pattern of thin wires. Given the mesh geometry, we must determine 
the magnitude of negative bias voltage necessary to maintain a repelling elec- 
tric field at the lunar surface. I set up a simple three-dimensional simulation 
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Fig. 5. — -Photoelectric yield Y for neutral graphite and silicate grains as 
a function of incident photon energy hv for several values of the grain size 
a, as indicated. 

Figure 5: Photoelectric yield as a function of photon energy for graphite and 
silicate grains 
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with the HiPhi code to provide the information. The simulation geometry, 
shown in Fig. 6, consists of an infinite mesh of thin wires in a square pattern 
with side length D. A ground plane is located 2 m below the mesh. A plane 
2.0 m above the mesh is biased to 5.0 MV to create a gradient E z = 2.5 
MV/m. Because of symmetry it was necessary only to model one quadrant 
of a mesh cell. The Neumann boundary condition (E parallel to the surface) 
was applied on the four sides in the x — y plane. The shielding effectiveness 
depends mainly on D relative to the distance to the ground plane - there is 
little dependence on the wire diameter. Figure 7 shows calculated equipoten- 
tial lines near the ground shield in the plane y = 0.0 with D = 40.0 cm and 
no bias (shield at ground potential). If the wires have diameter W — 0.32 
cm (1/8”), the opacity of the mesh is 


AWD/2 2 W 

r = ~d^ = ~d 


0.032. 


( 3 ) 


The field partially penetrates the mesh - the average potential beneath 
the wire array is about +64 kV. As a result there is an approximately uniform 
field E z = —33.4 kV/m in the region between the ground shield and lunar 
surface, sufficient to extract all photoelectrons. To ensure electron repulsion 
it would be necessary to apply a bias voltage of about -70 kV to the ground 
shield. Reducing the wire spacing to D = 20 cm doubles the opacity but 
significantly reduces the field in the protected region. With a grounded 
mesh, the field is E z = 6.3 kV/m, so a bias voltage of about. -15 kV would be 
sufficient. A double-layered ground shield would be desirable for personnel 
protection. The ideal ground shield system may consist of a grounded mesh 
at sufficient height to permit operations and a second mesh 1-2 m above with 
negative voltage to ensure E z has a positive value at the lower mesh. 


5 Shield field configuration and stored energy 

I set up a simulation to illustrate the field geometry of a finite-dimension 
shield system and to calculate the stored electrostatic field energy. I assumed 
cylindrical symmetry and applied the two-dimensional EStat code. Figure 8 
shows the simulation geometry. The moon surface at the left is treated as 
infinite ground plane. The ground shield 2 m above the surface has a 34 m 
radius and is biased to -60.0 kV. The HV shield 10 m above the surface is 
biased to +20 MV. The outer radius of the grading electrode is 21 m. I placed 
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Surface quantity: tz 

Plane position: Z = -4.000E-G2 
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Figure 6: Geometry - HiPhi simulation of field penetration of an infinite 
mesh with wires on a square patter with D = 40 cm. The colored plane 4.0 
cm below the wires shows E z . 
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Figure 7: Equipotential lines in the plane y = 0.0 cm near a mesh with 
D = 40 cm. Ground plane to the left, HV shield to the right. 
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Figure 8: Simulation of a finite-dimensional shield system - geometry and 
equipotential lines 

a grounded boundary at z = 100 m and r = 100 m to represent infinite space. 
Figure 8 shows calculated equipotential lines - electric fields are normal to 
the lines. Taking an integral of eoE 2 /2 over the solution volume gives a 
stored electrostatic energy of 0.68 MJ. This amount of energy could cause 
significant damage in the event of a breakdown, emphasizing the necessity of 
the ground shield. 

The field lines in Fig. 8 have a significant radial extent. Figure 9 shows 
a scan of E z in radius near the lunar surface (z = 1.0 m). The strong 
negative field values outside the radius of the ground shield could accelerate 
photoelectrons from a large area. Clearly, the radius of the ground shield 
must be substantially larger in a practical system. 
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Figure 9: Radial scan of E z close to the lunar surface (z = 1.0 r 
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6 Estimating shield efficiency 

I could perform extensive three-dimensional ion orbit calculations to estimate 
the effectiveness of electrostatic shields of the type of Fig. 2. This effort 
would be largely unnecessary because it is straightforward to estimate an 
upper limit on the shield efficiency for all configurations. Because the shield 
does not change the kinetic energy of penetrating ions, I define the efficiency 
E s as the ratio of the reflected ion flux to the incident flux. For solar proton 
events which could provide a lethal dose, the quantity E s must be close to 
unity. The following assumptions lead to simple model: 

• The incident ions have an isotropic distribution (equal flux per solid 
angle) . 

• Ions arrive from free space at the same reference potential as the lunar 
surface. 

• The shield is a plane of infinite transverse extend biased to +V 0 - 

Regarding the third condition, finite assemblies like Figs. 1 and 2 have regions 
with <p < Vo and therefore have lower efficiency than the ideal system. 

The polar angle 9 is defined with respect to the z axis in Fig. 8. It equals 
0.0° for normal incidence and 90.0° for ions arriving from the horizon. With 
the assumption of equal flux per solid angle, the normalized probability as a 
function of angle is given by: 


p(9)d9 = sin (9)d9. (4) 

Consider a non-relativistic ion incident from oo at angle 9. The ion has kinetic 
energy T, ionization state Zi and kinetic energy per unit charge T 0 = T/Zi. 
The component of kinetic energy normal to the lunar surface is 


T± = T 0 cos9. (5) 

Because the electric fields are normal to the surface, the condition that 
the shield reflects an ion is that eVo > T±, or cos# < cos#o, where: 


cos#o = 



(6) 


32 



4 



Figure 10: Upper limit on shield efficiency as a function of ion kinetic energy 
per charge for Vo = 20 MV 

Integrating Eq. 4 from 9 — 90.0° to 0.0° gives the fraction of reflected ions, 
equivalent to the shield efficiency: 

I eV n 

E s - sin 6d6 = cos 9 0 = \ — . (7) 

Jn/2 V T 0 

To illustrate the implications of Eq. (), suppose the shield bias is +20 MV. In 
this case, the shield reflects all ions with kinetic energy per change To < 20.0 
MeV. At T 0 = 40.0 MeV, the reflection coefficient is E s = 71%. Figure 10 
shows E s for the choice Vo = 20.0 MV over the energy range from To = 10.0 
MeV to 1000.0 MeV. Clearly, the shield is ineffective for biological protection 
from the high-energy flux of gradual solar proton events (~ 100 MeV) and 
galactic cosmic radiation (~ 1 GeV). 
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